Long-Range Charge Transfer in DNA through Polaron Diffusion 
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Recent experimental evidence shows that the ir orbitals along the stacking of base pairs can 
facilitate the long-range charge transfer in DNA [5-7]. Proton motion in the base pair hydrogen 
bonds has also been found to affect the transfer rate. To explain this behavior we propose a model 
considering interactions of doped charges with hydrogen bonds and vibrations in DNA. The charge 
trapped by either protons or vibrations can cause structural distortions leading to polaron formation. 
By further considering polaron diffusion in DNA we find that the charge transfer rate derived from 
the diffusion coefficient is in agreement with the experimental results [33]. 



I. INTRODUCTION 

The question on whether DNA is a conductor or an 
insulator is still controversial due to the different results 
obtained in the measurements of DNA conductivity . 
The related issue of long-range charge transfer in DNA 
which is associated with the problem of charge mobil- 
ity has also been discussed for many years with an equal 
amount of heated debate. It has been suggested that 
the overlap of the electronic ir orbitals along the stacked 
base pairs provides a pathway for charge propagation 
over 50 A and the reaction rate between electron donors 
and acceptors does not decay exponentially with the 
distance 5-7 . 

A multiple-step hopping mechanism was recently pro- 
posed to explain the charge transfer behavior in DNA 8 . 
In this theory, the single G-C base pair is considered as 
a hole donor due to its low ionization potential if com- 
pared to the one on A-T base pairs. Long-range charge 
transfer is accomplished by a series of incoherent charge 
tunneling events between two nearest G-C base pairs sep- 
arated by A-T pairs. The probability of the short-range 
tunneling (super-exchange) strongly depends on the dis- 
tance between the two G-C base pairs. However, recent 
experimental results for the reaction rate only showed a 
very slight distance dependence for the charge hole trans- 
fer through more than three A-T base pairs 9 . It is now 
believed that a hole can be created in a A-T base pair by 
thermal activation from a G-C base pair when charge 
tunneling becomes unlikely over a long distance 10-12 . 
The long-range hopping can, therefore, be carried out 
through a long A-T bridge just as by a series of short- 
range tunnelings in the G-C pairs. Although this mech- 
anism seems to give a good interpretation of the charge 
transfer in DNA, the cause of the incoherent charge hop- 
ping over localized states of the hole carriers with such 
low reaction rate (10 9 ~ 10 6 s -1 ) has not yet been well 
understood. 

Besides the multiple-step hopping mechanism, polaron 
motion has also been considered as a possible mechanism 



to explain the phenomena of DNA charge transfer 13 ' 1 . 
The charge coupling with the DNA structural deforma- 
tions can create a polaron and cause a localized state. 
From the study of the dynamical properties in one dimen- 
sional systems, it is known that the polaron behaves as 
a Brownian particle that collides with the low energy ex- 
citations of its environment which acts as a heat bath 15 . 
This diffusive behavior is very similar to the multiple- 
step hopping mechanism based on a random walk model, 
although the former occurs in a continuous media while 
the latter is considered on the discrete sites of a lattice. 
Thus, the behavior of the incoherent charge hopping can 
be understood as polaron diffusion. 

For most physical systems acoustical and optical 
phonons are the main lattice excitations. However, when 
a system contains hydrogen bonds, proton motions also 
need to be considered. Instead of oscillatory motions in 
a single-minimum potential, protons can tunnel from one 
side of a hydrogen bond to another in a double-minimum 
potential 16 . This proton tunneling causes interstrand 
charge hopping and is important in spontaneous muta- 
tions. Since genetic information can be preserved only 
when a nucleotide base is matched with its complemen- 
tary one by hydrogen bonding, the positions of protons 
are the main identification for DNA polymerases to make 
replications with high fidelity. However, the reaction of 
proton transfer could generate tautomeric base pairs and 
destroy the fidelity 17 . The mechanism of the proton tun- 
neling in the isolated base pairs and its possible biological 
implications were widely discussed by Lowdin many years 
ago 18 . It has also been suggested that proton transfer and 
charge conduction in stacked base pairs are affected by 
each other 19-21 . The recent experimental results on the 
influence of mismatched base pairs show that the proton 
transfer is required for the long-range charge transfer 22 . 
Thus the effect of proton tunneling will be discussed here. 

The paper is organized as follows. The model for 
charge motion in DNA will be introduced in the next Sec- 
tion, considering charge hopping, possible fluctuations 
in structure and their interactions. We will present the 
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formation of a polaron in Section III, together with the 
study of polaron stability and proton delocalization in the 
hydrogen bonds. In Section IV, we calculate the polaron 
diffusion coefficient and the results will be compared with 
the experimental reaction rate of charge transfer in DNA. 
Section V contains our conclusions. 



II. THEORETICAL MODEL 

Although biologically functional DNA is an aperiodic 
system, we will focus here on a periodic DNA structure. 
Since the polaron state discussed here is highly localized 
(see below), we believed that the periodicity condition is 
not very constraining when we consider charge transfer 
over intermediate distance. In fact, in some experimen- 
tal setups one looks for charge transfer between two G-C 
pairs which are separated by many (periodic) A-T pairs. 
As we can see in Fig. 1, following hole injection into the 
system, the charge localizes in the G-C base pair but can 
migrate from one unit cell to another. The essential in- 
gredients in our model are (i) the electronic system, that 
realizes the charge conduction; (ii) the possible structural 
fluctuations, which include phonons and proton motion 
in hydrogen bonds; and (Hi) the coupling between the 
charges and the structure. 
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FIG. 1. Charge hopping along the stacked base pairs of a 
periodic DNA chain, m is the number of the charges in the 
i-th unit cell. 

For a non-interacting electron system, the tight- 
binding model is extensively used in many theories of 
DNA charge transport 11 ' 14 ' 23 : 



H e = £ -t (C}C i+1 + C\ +l C t ) 



(1) 



4io , where a is the lattice spacing and the wave vectors 
k is 2im/L with n as an integer and the system size L. In 
the continuum limit (a — > 0), the Hamiltonian describes 
the kinetic energy of charge particles with a Bloch wave- 
function and effective mass m , where h 2 /2ma 2 is used 
to replace the overlap integral to of the discrete model. 

For the phonon part of the structural fluctuations in 
DNA, we consider that the relative motions between two 
different base pairs can be represented by an acoustical 
phonon mode and the vibrational motion inside a base 
pair by optical phonons 25 . These modes represent lattice 
distortions, such as sliding, twisting or bending. The 
Hamiltonian that describes these modes is: 



H p h — 
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2M + i Mt ^ ( Ui+1 ~ Ui ^ 



(2) 



where Ui and Vi are the lattice displacement and the in- 
ternal vibration coordinate of the i-th unit cell, pi and Pi 
are their conjugated momentum, respectively, M is the 
mass of a unit cell, wo the oscillation frequency of the op- 
tical phonon and the dispersion relation of the acoustical 
motion is: lo s (k) = c s k , where c s is the sound velocity 
along the chain. 

The charge coupling to acoustical phonons is given 
by the Su-Schricffcr-Hecger model 2 and the interaction 
with optical phonons is described via the molecular crys- 
tal model of Holstein 26 . Thus the total electron-phonon 
interactions of DNA system is 



(3) 



where j s and j v are the coupling constants. In sys- 
tems with half-filled conduction band, such as an one- 
dimensional 7r-conjugated polymers, the SSH term gener- 
ates dimerization in the ground state (Peierls instability) 
and forms solitons in the excited states 27 . However, for 
DNA, considered here as a band insulator, both interac- 
tions can generate lattice distortions and lead to polaron 
formation when a charge is doped into the molecule. The 
second term in (3) usually induces small polarons in ionic 
crystals. 

For the proton motion in the hydrogen bonds we use a 
two-level system to describe the tunneling behavior (see 
Fig. 2) 28 : ' 



where Ci(Cj) is the electron annihilation(creation) oper- 
ator at the i-th site and to is the overlap integral of elec- 
tronic 7T orbitals between two nearest neighboring sites. 
(We disregard the electron spin since we are not describ- 
ing magnetic phenomena here.) The energy spectrum of 
the electrons is Tiujk = — 2t cos (ka) giving a bandwidth 



H a = -(-ea z +ta x ). 



(4) 



where a z and a x are Pauli matrices, e is the energy bias 
between the two localized proton states and t is the tun- 
neling matrix element. The ratio t/e provides informa- 
tion of the proton motion in the system. When the ratio 



2 



is small, the protons are localized in one side of the hy- 
drogen bonds. However, when it is large, the protons 
are delocalized. As we can see in Fig. 2(a), the nor- 
mal G-C base pair is the lower energy state and its tau- 
tomeric form (G*=C*) is an excited state. When e 3> t , 
the probability of having a tautomeric form is extremely 
small. However, in Fig. 2(b), the radical cation of a G-C 
base pair has almost the same energy as its tautomeric 
form 29 , i.e., s < t . In this case, the proton state becomes 
delocalized in the hydrogen bonds. 



around the hydrogen bond (n, = 0), the coupling van- 
ishes. Although there are many hydrogen bonds in the 
base pairs of a unit cell, it is not necessary to consider 
all of the possible proton transfers since at room tem- 
perature the transitions to the tautomeric forms are very 
unlikely to occur. We only need to consider the hydrogen 
bond with the most probable proton transfer with charge 
injection into the unit cell, i.e., only one two- level system 
is required for each unit cell. 

By giving a total Hamiltonian as the sum of (1), (2), 
(3), (4) and (5) we have our complete model for charge 
motion in DNA. 
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FIG. 2. Proton transfer in the hydrogen bond of a DNA 
base pair described by a two-level system, (a) For a neutral 
G-C base pair, proton transfer is unlikely to cause the tau- 
tomeric form and we have £ >i . (b) The radical cation of 
a G-C base pair only has a little energy difference with the 
tautomeric form, i.e., e < t. 

To model the coupling between the protons in the hy- 
drogen bond and the charges in the DNA strand, we use: 



H e -a = la (of - l)rn 



(5) 



where 7 CT is the coupling intensity and rii — CjCi is the 
number of charges on the i-th site. When the proton is 
in the lower energy state (of = 1) or there is no charge 



III. FORMATION OF POLARON IN DNA 

Although polaron formation has already been consid- 
ered very important for charge motion in DNA 13,14 ' 30 , 
the effect of the proton motion in the hydrogen bonds has 
not yet been discussed in the literature. Here, to investi- 
gate this factor, we analytically solve our model by tak- 
ing the Hamiltonians into the continuum limit. The car- 
rier wavefunction ip (x, t) and the phonon modes u s (x, t), 
u v (x, t) are considered as functions of time and a contin- 
uous coordinate x .The equations of motion for the vari- 
ables in the problem can be studied in the semiclassical 
limit (h — > 0) and by a change of variables A = x — vt , 
where v is the polaron velocity. In this case, we have 15 : 



u s (v) (x, t) = u,(„) (a; - vt) = u s ^ (A) 
tp (x, t) = 0o (A) exp [j- (mvx - E t)] , 



(6) 



where E is the binding energy of the polaron, and the 
problem reduces to a nonlinear Schrodinger equation for 
the wavefunction </>o : 



h z d 2 



2m d\ 2 



E 



j a sin" 8 



+ 



p{c 2 -v 2 ) puj 2 ) 







(7) 



where p is the mass density M/a , sin 2 6 represents the 
positions of the protons in the hydrogen bonds and is 
expressed by 



sin 2 9 



1 



1 



Vs^+l 2 



(8) 



where e = £—jo-(f>o 2 ■ The interaction energy between the 
carriers and protons, j a sin 2 9 , is significant only when 
the charge trapped by the hydrogen bond causes proton 
derealization. Therefore, for small value of t compared 
to e , we can use the approximation: 
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i(e-t) 



-t < 7<r 
, 7 CT ^o 2 > e + t 



<e-t 
fo 2 <£+t 



(9) 
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From the continuity of the polaron wavefunction, equa- 
tion (7) can be analytically solved and the solution is 
found as combinations of Jacobian elliptic functions. 

When the charge coupling with the protons is small, 
i.e., 7 CT (</>o)Lm — £ ~ t (where (0o) max is the maximum 
value of 4>o), the polaron is purely formed by the inter- 
action with phonons and the proton motion has no effect 
on the wavefunction or the binding energy The ground 
state solution is simply expressed as 



<t>Q 0) = V9 sech (2g\) , 



(10) 



and the binding energy is 



pair in the center of the polaron, its tautomeric form can 
be easily induced. 

The relation between the coupling constant 7 CT and the 
binding energy is shown in Fig. 3, where we have 1.7 eV 

i 2 2 

A for the value of the coupling strength + . When 

j a is larger than 0.24 eV , the proton becomes delocalized 
in the hydrogen bond and the charge binding energy in- 
creases rapidly. The local deformations from the charge 
interaction with phonons can also stabilize the polaron 
as shown in Fig. 4 with 7 CT = 0.25 eV. The increase of 
the couplings 7 S or j v also induces proton transfer since 
the carrier is strongly attracted in the base pair. 
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where 



m 



+ 



p(c 2 8 -v 2 ) pul 



(11) 



(12) 



-E g (tV) 



It is easy to see that l/g represents the size of the po- 
laron. To calculate the binding energy and the size of a 
static polaron, we consider that a hole carrier is doped 
into a poly(G)-poly(C) DNA. The effective mass of the 
carrier™ = fi /2t Q a 2 is of order 20m e = 1.8219 x 10~ 29 kg 
(m e the mass of a electron) obtained from band structure 
calculations, which give a bandwidth of about 0.04 eV 31 . 
For the charge coupling strength ^» _|_ 1^ generated 

by the acoustical and optical phonons in (7), since the 
energy difference between the adiabatic and vertical ion- 
ization potentials of a G-C base pair is about 0.5 eV 29 , 
the optical part is given as about 0.5 eV x3.38A = 1.7 
eV A. The acoustical part is assumed to be smaller, so 
we consider the value of the coupling strength as approx- 
imately 2 eV A. From these results, we find that the 
polaron is 0.75A wide with binding energy E = —1.35 
eV and, therefore, it is entirely localized in a base pair. 
For the weak coupling to hydrogen bonds (7^ < 0.2 eV), 
the protons are localized in their original positions and 
the tautomeric base pair is unlikely to be generated. 

As the coupling constant 7 CT becomes larger, proton 
transfer can be induced since the charge trapped into 
a base pair reduces the bias energy e dramatically. By 
considering the bias energy e and the tunneling matrix 
element t of order 0.3 eV A -1 and 0.03 eV A -1 , respec- 
tively, we find that the binding energy of the polaron Eg 
is —1.45 eV for -f a = 0.226 eV, which is more stable than 

the above with the same values of the effective mass m 

2 2 

and the coupling strength ^ + . The ground state 

wavefunction consists of two parts, the central and the 
external. For the external part, the wavefunction behaves 
as a hyperbolic secant function since it should vanish as 
A — > 00 . In the central part of the polaron, the wave- 
function is described by a Jacobian elliptic functions of 
second kind and the proton in the hydrogen bond is de- 
localized since the value of sin 2 9 is large. For the base 
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FIG. 3. Polaron binding energy as a function of the cou- 
pling constant between charges and protons 7 CT , where the 
parameters are given by: ^ = 2.624 eV _1 A -2 , e = 0.3eV 
A -1 , t = 0.03 eV A~ x and the charge coupling strength with 
phonons is 1.7 eV A. 
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FIG. 4. Polaron binding energy as a function of the cou- 

2 2 

pling strength between charges and phonons — P% — ^- + 7 " ■> 

Ps(c 2 -i> 2 ) Pv"£ 

, where we use the same parameters as in Fig. 3 and the 
charge coupling constant with proton 7 CT is 0.25 eV. 

The effect of proton transfer has been observed in ex- 
periments on the influence of mismatches in the DNA 
long-range charge transfer 22 . Mismatched base pairs are 
found to cause a large drop in the charge transfer rate due 
to the weakening of hydrogen bonds. For a normal G-C 
base pair, the positive charge in the radical cation can be 
distributed by both nucleotides through the delocaliza- 
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tion of protons. However, mismatched Guanine, which 
leads to a fast deprotonation to the surrounding water, 
might induce a bound state with lower energy. There- 
fore, the charge is localized in the mismatched base pair 
and the long-range charge transfer is interrupted. 



IV. POLARON DIFFUSION IN DNA 

The dynamic properties of a polaron induced by the 
carrier interaction with acoustical and optical phonons 
in one dimension were explored many years ago 15 . In the 
strong coupling limit, it is found that the polaron moves 
as a Brownian particle interacting with the light parti- 
cles of the environmental heat bath. To investigate these 
results in the context of DNA, we first consider the case 
without proton coupling (i.e., 7 CT = 0). 

The damping parameter of the optical polaron motion 
at high temperature (fcsT ^> %u)q) is given by 15 : 



= 7T k B T g a 2 
lo 6 oj q M>' 

where M is the classical polaron mass: 



M„ = 



rn 



Eg 



and 



9o = 



m jy 
h 2 puj 2 



~~ ft 5 " " 



(13) 



(14) 



(15) 



Here, we use E as the binding energy from the optical 
phonons. For acoustical phonons we have the damping 
parameter: 



la 



315 k B T g a 



32tt 4 



M ' 



(16) 



where the classical polaron mass for the acoustical case 
is 



M n = 



32m / E, 



frc s g a 



and 



g a 



m 7 S * 
h 2 p s c 2 



(17) 



(18) 



with E a the binding energy from the acoustical phonons. 
The random walk motion of the polaron in the long time 
limit (t ^> l/7 or t ^> l/7 a ) is governed by the one- 
dimensional diffusion equation: 



dn £)d 2n 
dt dx 2 ' 



(19) 



where n (x, t) is the density of the particles and D is the 
diffusion coefficient. The relation between the damping 



parameter and the diffusion coefficient is given by the 
Einstein relation: 



D (a) = 



lo(a) M o(a) 



(20) 



where we use -D ( a ) for the diffusion coefficient in the 
optical (acoustical) case. 

As we can see from the scheme of the multiple-step 
hopping mechanism in Fig. 1, the hole hopping from a 
G-C base pair to its nearest neighbor can be viewed as 
an oxidation-reduction reaction with a rate constant k . 
According to the rate law of a first-order reaction, the 
concentration of the charge at the i-th G-C base pair, n i} 
is described by 



drii 
~dt 



= k (n i+1 + Hi-i - 2rii) . 



(21) 



In the continuum limit, when the distance between two 
nearest neighboring G-C base pairs, a, is very small, we 
have 



dm 
~dt 



= ka 2 



(n i+1 -ni) 



ka 



2 d 2 n l 
dx 2 ' 



(22) 



By comparing this equation with (19), we have the corre- 
spondence from the diffusion coefficient in the continuous 
medium to the charge transfer rate of the discrete model: 



D 



ka 2 . 



(23) 



Thus, instead of a diffusion coefficient, we can use the 
discrete characteristic to to estimate the rate constant. 
From (20) and (13), we have 



D n = 



k B T 6 Wq 3 uiq fi 



J M 7r g a 2 2ir E Q 2ma 2 



6 0J 2 

^E~ o t0a ■ 



(24) 



By the correspondence (23), the rate constant for the 
optical case is obtained as 



i, 3 *° 
fco= 2^°- 



(25) 



We can also derive the rate constant for the thermal ac- 
tivation by the acoustical phonons as 



64tt 4 Ps c s 3 t 
315 7s 2 ' 



(26) 



For poly(G)-poly(C) DNA, the oscillation frequency 
is ujq w 10 11 Hz from the theoretical and experimental 
results 32 . The theoretical calculation has shown the va- 
lence bandwidth of 0.04 eV 31 . By considering that the 
binding energy E Q is about 1 eV, we have the reaction 
rate k Q ~ 10 9 s _1 . This estimate is in agreement with the 
experimental result for electron transfer in DNA which 
gives k a w 10 8 - 10 9 S - 133 . 
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This low reaction rate is generated by two factors in 
equation (25), the hopping probability J?- and the oscilla- 
tion frequency loq . Due to the strong structural fluctua- 
tions at room temperature, the charge hopping could only 
happen when the local structure of the polaron achieves 
an effective configuration 30 . For each period of structural 
oscillation, the probability to have a successful hopping is 
of order J?- . Thus, the charge transfer rate is determined 
by the frequency of reaching the effective configuration 
iv and the probability The behavior of the inco- 
herent charge hopping is now described by the polaron 
motion under thermal fluctuations. 

Another way to study the diffusion cocfficcnt in the 
optical case is to consider the mean square displacement: 

. poo 2 

(Ax) 2 ) = (x 2 ) = / x 2 <f> (x) O (x) dx = — , (27) 



and we may rewrite (24) as: 



D n 



72 lu 



(AxY 



(Ax)' 



TT 2 27T 



2tt 



(28) 



Now, we can see that, for every oscillation period 2n/u>o 
, the collisions by the light particles from the environ- 
ment spread the polaron with the mean-square displace- 
ment (^(Ax) 2 ^. The diffusion in the continuous medium, 

however, needs to be corrected when we make the cor- 
respondence since the size of the polaron is smaller than 
the lattice spacing. Therefore, the real reaction rate of 
charge transfer between the discrete base pairs should 
be smaller than our results. Nevertheless, the order of 
the charge transfer rate can still be estimated from the 
simple expression: 



k 



w 



(Ax)" 



2tt 2a 2 



(29) 



We can also estimate the reaction rate k a in the acous- 
tical case from the experimental data for the sound veloc- 
ity c s « 10 5 cm s 32 . Assuming ~ -34 eV A, we have 

k a « 10 13 s -1 which is four order of magnitude larger 
than k . A faster charge transfer rate is expected if the 
polaron is purely generated by the acoustical phonons. 
However, if optical phonons or proton transfer arc in- 
volved, a smaller value of the reaction rate should be 
expected clue to the increase of the binding energy in 
those cases. 



V. CONCLUSIONS 

In this paper we studied the charge conduction in pe- 
riodic DNA focusing on the effect of the structural fluc- 
tuations which includes the proton motion in hydrogen 
bonds. We found that the charge trapped by hydrogen 
bonds or phonons can generate a polaron state and the 



coupling to the protons can induce proton transfer in a 
base pair. The polaron moving with small velocity is seen 
as a Brownian particle colliding with the light particles of 
the environment. This diffusion process corresponds to a 
multiple-step hopping mechanism in the discrete model. 
From this correspondence, the reaction rate of the long- 
range charge transfer in DNA can be derived from the 
diffusion coefficient, and it was found that the result pre- 
dicted in the optical case is in agreement with the exper- 
imental results. 

Although we did not include the effect of proton trans- 
fer on the polaron diffusion, it is expected to be crucial in 
the dynamics of charge hopping. From electrochemical 
experiments 34 , a decrease of charge transport along the 
DNA with mismatched base pairs is found. This behavior 
might be explained by the effect of the proton motions in 
the hydrogen bonds and we will explore this possibility 
in future work. 
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